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摘 要 : 巴 丹 吉林 沙漠 腹地 分 布 有 大 量 湖泊 ,水 面 薰 发 是 湖泊 水 循环 万 至 整个 沙漠 水 循环 的 关键 控制 因素 。 本 文 
以 在 沙漠 腹地 湖泊 一 一 苏 木 巴 润 吉林 建立 自动 气象 站 监测 的 2012 年 9 月 11 日 至 2013 年 3 月 26 日 和 2013 年 5 H 


8 日 至 2015 4E 6 H 23 日 的 长 序列 气象 落 发 数据 ,以 及 2012 年 9 月 13 日 至 2014 年 6 
蒸发 与 气象 要 素 的 动态 关系 。 利 用 多 元 线性 回归 模型 非 线性 空气 动力 学 模 
蒸发 量 进行 模拟 分 析 。 结 果 表 明 : 随 着 时 间 尺 度 的 提升 , 落 发 量 与 同期 气 


基础 ,在 多 时 间 尺 度 下 分 析 湖 泊 水 耳 
型 和 Penman 公式 对 逐 小 时 、 逐 日 和 逐 


120 日 的 湖底 压强 数据 为 


象 要 素 之 间 的 相关 性 增强 。 在 逐 小 时 尺度 下 ,蒸发 与 同期 气象 要 素 的 相关 性 较 弱 ,用 多 元 线性 回归 模型 拟 合 效 
不 佳 。 在 逐日 和 逐 月 时 间 尺 度 下 ,蒸发 量 与 同期 气温 和 水 温 的 相关 性 最 高 , 利 


FH, 
ps 
E 
气 


j 多 元 线性 回归 模型 或 非 线性 空 


动力 学 模型 可 以 较 好 地 模拟 蒸发 量 。Penman 公式 计算 的 逐日 敬 发 量 在 夏季 偏 大 而 冬季 偏 小 。 人 研究 结果 可 为 沙漠 


湖泊 水 量 蒸发 的 评估 提供 科学 依据 。 
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水 文 循环 的 研究 需要 对 气象 要 素 (如 降水 L2 
发 等 ) 和 水 文 要 素 ( 如 径流 、 人 人 渗 等 ) 进行 直接 观测 
或 间接 评估 ,前 人 已 经 发 展 了 很 多 量化 方法 "5 。 
这 些 要 素 中 ,降水 和 径流 的 观测 技术 相对 成 熟 可 靠 、 
而 蒸发 的 准确 测量 则 较为 困难 !9 DL EC ZR f t E 
化 成 为 水 文学 研究 的 难点 和 热点 "”] 。 准 确 计 算 水 
体 莹 发 对 湖泊 和 湿地 水 资源 及 生态 系统 的 评估 管理 
至 关 重 要 ,尤其 是 在 水 资源 比较 稀缺 的 地 区 中。 

巴 丹 吉 林 沙 漠 位 于 我 国 西北 干旱 地 区 ,沙漠 腹 
地 湖泊 众多 ,其 中 常年 有 水 的 湖泊 有 110 48027, 
由 于 降水 量 与 水 面 莹 发 量 的 巨大 收 支 差 距 ,湖泊 的 
形成 机 理 和 演变 趋势 成 为 重点 关注 的 科学 问 
WUO 。 最 近 十 几 年 来 ,大 量 学 者 在 巴 丹 吉 林 沙 漠 
进行 了 科学 探究 调查 ,但 是 对 于 这 一 问题 仍然 没有 
达成 共识 。 由 于 该 沙漠 没有 地 表 径 流 , 蔡 发 几乎 是 
湖泊 唯一 的 排泄 方式 ,也 是 沙漠 水 循环 的 关键 环节 。 
因此 要 了 解 湖泊 的 形成 机 理 , 就 必须 先 研究 清楚 湖 
泊 的 落 发 强度 及 其 如 何 变 化 。 前 人 对 湖泊 的 蒸发 量 
进行 了 估算 ,但 不 同 结果 之 间 存 在 较 大 差异 。2004 
年 ,Chen 等 "7 认为 巴 丹 吉 林 沙 漠 腹 地 的 湖泊 年 蒸 
发 量 可 达 4 000 mm; 2008 年 , Gates 等 "推测 该 地 
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区 湖泊 的 潜在 蒸发 量 为 2 600 mm * a^ ' ; 2012 年 ， 
马 宁 ' ”使 用 涡 动 相关 技术 进行 观测 ,根据 能 量 平 衡 
原理 计算 出 湖泊 的 年 燕 发 量 约 为 1 400 mm;2014 
年 , 王 旭 升 等 "在 巴 丹 吉 林 沙 漠 进 行 了 水 文 地 质 调 
查 和 典型 湖泊 综合 观测 ,分 析 认 为 沙漠 湖水 的 蒸发 
强度 为 1 200 ~1 500 mm * a^! , HE SF. Yang 等 
基于 Penman 公式 计算 得 出 的 沙漠 腹地 湖泊 年 蒸发 
ft (1040 mm) ;2015 年 , 胡 文 峰 等 .1 利用 2012 年 8 
月 至 2013 年 7 月 在 巴 丹 吉林 沙漠 腹地 音 德尔 图 湖 
泊 上 用 涡 动 相关 系统 观测 的 潮流 资料 ,计算 得 到 湖 
泊 的 年 燕 发 量 为 (1 450 + 10) mm;2017 4E, 5K 3e 
等 引 通 过 在 苏 木 吉 林 湖 区 建立 地 下 水 - 湖水 动态 
模拟 模型 , 拟 合 2012 一 2013 年 的 湖水 位 季节 性 波 
zh ,推测 湖面 的 年 蒸发 量 为 1 261 mm。 这 些 研究 主 
要 局 限 在 年 尺度 的 湖面 蒸发 量 评估 ,对 日 尺度 和 月 
尺度 的 蒸发 特征 分 析 较 少 。 

从 机 理 上 讲 , 巴 丹 吉 林 沙 漠 湖 泊 的 蒸发 量 及 其 
动态 特征 受到 太阳 辐射 ,空气 动力 学 条 件 和 湖水 传 
热 等 因素 的 影响 。 在 不 同时 间 尺 度 上 ,哪些 因素 对 
湖面 蒸发 起 主导 作用 ,如 何 根据 气象 要 素 的 观测 推 
算 湖面 荧 发 强度 ,这 是 有 竺 解决 的 实际 问题 。 本 文 
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利用 在 巴 丹 吉林 沙漠 腹地 典型 湖泊 进行 监测 得 到 的 
长 时 间 序 列 气 象 数据 和 落 发 上 数据 ,分 析 燕 发 与 气 
象 要 素 的 动态 关系 , 旨 在 获得 不 同时 间 太 度 下 湖泊 
蒸发 量 对 气象 条 件 的 依赖 性 ,为 研究 巴 丹 吉 林 沙 漠 
水 循环 提供 新 的 参考 。 


1 研究 区 概况 


巴 丹 吉林 沙漠 (39"30' ~ 4220' N,99?30' ~ 
104*0'E) 属于 内 蒙古 阿拉 善 盟 ,位 于 黑河 流域 下 游 
盆地 的 东部 , 北 临 拐子 湖 , 东 接 宗 乃 山 ,东南 和 西南 
分 别 以 雅 布 赖 山 和 北大 山 为 界 ,总 面积 5.2 x 10* 
kmh ,属于 中 国 第 二 大 沙漠 (图 1) 。 该 沙漠 的 地 
势 总 体 上 东南 高 西北 低 ,沙漠 腹地 是 高 大 密集 的 沙 
丘 和 沙 山 , 沙 山 面积 可 占 沙 漠 面 积 的 60% 以 上 ,其 
相对 高 度 一 般 在 200 ~ 300 m。 巴 丹 吉 林 沙 漠 处 于 
亚洲 季风 西北 边缘 区 , 主要 受到 中 纬度 西风 带 和 亚 
洲 季风 的 双重 影响 , 冬 春季 主要 盛行 西风 西北 风 ， 
夏秋 季 主 要 盛行 东南 风 , 气 候 干 旱 ,多 年 平均 降水 量 
呈现 东南 多 西北 少 的 规律 ,潜在 蒸 散 量 从 南 到 北 增 
加 。 地 下 水 的 总 体 流 向 是 由 南 向 北 .由 东 向 西 ,这 与 
该 沙漠 的 总 体 地 势 南 高 北 低 、 东 高 西 低 相 一 致 ,其 水 
力 梯度 为 0.8%o ~ 7.9%o05] 。 巴 丹 吉 林 沙 漠 腹 地 湖 
泊 众 多 ,但 湖泊 面积 一 般 都 小 于 0.2 km ,其 中 最 大 
的 是 诺尔 图 ,面积 达到 1.5 km , 除 诺尔 图 外 ,还 有 4 
个 湖泊 面积 大 于 1 km ,包括 苏 木 巴 润 吉林 、 呼 和 吉 
林 、 音 德尔 图 与 车 日 格 勒 。 沙 漠 湖 泊 绝 大 部 分 为 咸 
水 湖 , 矿 化 度 在 1~400 g - L”。 本 文 观测 研究 的 湖 
泊 是 第 2 大 盐湖 , 即 苏 木 巴 润 吉林 , 面积 》 
1.24 kn? ,最 大 深度 超过 11 mO (图 1)。 


2 监测 数据 


为 了 监测 苏 木 巴 润 吉林 的 气象 要 素 , 项 目 组 在 
湖 中 搭建 了 钢 架 木板 桥 平 台 ,用 以 安装 自动 气象 站 
(图 1)。 该 气象 站 的 位 置 为 102.42°E,39.79°N, 主 
要 用 于 实时 监测 气温 (7 )、 水 温 (7s)、 相 对 湿度 
(Rj) KEE (Ap) G8 (W) 、 风 向 (W)、 净 辐射 
(RO KEP) 及 燕 发 租 内 的 水 位 。 支 架 顶 端 固 
定 有 风速 仪 .风向 标 及 净 辐 射 仪 ,风速 仪 . 风 向 标 距 
桥 面 高 度 为 202 cm, 净 辐射 仪 距 桥 面 185 cm, TERR 
桥 面 高 140 cm 处 安装 了 湿度 感应 器 。 其 下 是 数据 
采集 系统 。 在 水 平 距 离 支 架 135 cm 处 的 湖 中 ,使 用 
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1 研究 区 及 气象 站 位 置 


Fig. 1 Geographical location of the study area and the 


distribution of meteorological stations 


倒 锥 形 钢 架 固定 了 一 个 E601 79 Z& Az M, zi& IA 
的 水 取 自 浅 层 湖水 ,底部 安装 了 测量 精度 为 二 0.1 
mm 水 压力 传感器 ,用 于 自动 监测 水 位 的 变化 ,由 此 
可 以 计算 出 蒸发 量 。 每 隔 2 月 (夏秋 季 ) 或 4 月 ( 冬 
春季 ) ,对 营 发 亚 进 行 一 次 彻底 换 水 ,以免 盐 度 效 
nU 过 于 显著 而 影响 观测 结果 。 桥 下 安装 了 温度 
传感器 ,用 于 监测 湖面 以 下 10 cm .20 cm 处 的 水 温 ， 
由 数据 采集 系统 记录 数据 。 雨 量 计 安装 在 支架 附 
近 , 其 顶部 距 桥 面 73 em ,长 度 22 cm, 外侧 直 径 20 
cm。 在 靠近 湖底 的 固定 位 置 ,安装 了 一 个 CTD-Div- 
er 传感器 ,监测 湖底 附近 水 的 电导 率 、 温 度 和 压强 ， 
每 隔 2 h 自动 读 取 并 记录 1 次 数据 。 根 据 CTD-Div- 
er 记录 到 的 压强 变化 ,进行 气压 校正 ,可 以 计算 出 
iK CH) 的 变化 。 该 位 置 湖水 很 浅 , 深度 小 于 
0.6 m, 

本 文 收集 了 2 个 时 期 的 气象 观测 数据 ,第 1 个 
时 期 为 2012 年 9 月 11 日 13:30 至 2013 年 3 月 26 
日 10 :00 ,第 2 个 时 期 为 2013 年 5 月 8 日 17:00 至 
2015 年 6 月 23 日 21:00, 数 据 采样 间隔 为 30 min, 
湖底 压强 的 连续 数据 采集 期 为 2012 年 9 月 13 H 
0:00 至 2014 年 6 月 20 日 0:00, 数 据 采样 间隔 为 
2h, 图 2 以 2013 年 1 月 15 日 和 2013 年 8 月 15 日 
的 每 30 min 气温 、10 cm 处 水 温 (7; i0) 、20 em 处 水 
(Tso) ,相对 湿度 、 风 速 . 气 压 \ 净 辐射 , 蒸发 亚 测 
得 累积 蒸发 量 (Eo ) 和 每 2 h 湖底 压强 的 数据 为 例 ， 
来 显示 气象 数据 小 时 尺度 的 变化 特征 。 这 两 天 均 没 
有 降水 发 生 。2013 年 8 月 15 日 (夏季 ) 的 主要 特征 
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图 2 研究 区 2013 年 1 月 15 日 和 8 月 15 日 气象 观测 数据 
Fig.2 Meteorological data observed in the study area on January15 and August 15, 2013 


为 :平均 气温 为 26.59 C ,10 em 深度 处 的 湖水 平均 
温度 为 26.13 % ,与 平均 气温 相近 ,而 20 cm 处 平均 
水 温 则 为 24.83 C , 略 低 于 平均 气温 ;气温 比 水 温 的 
波动 强烈 ;平均 湿度 平均 净 辐 射 和 平均 风速 分 别 为 
33.17% .107.92 W - m ”了 和 3.26m:s ;8:30 空 
气 湿度 最 高 ,14:30 最 低 ,14:00 净 辐 射 最 大 ,23 :30 
最 小 ,二 者 一 天 内 波动 较 大 ;平均 气压 .平均 湖底 压 
强 分 别 为 873. 30 hPa 和 909. 84 hPa ; Z 7z IIIL*4 FH] 8 
IZEN 8.3 mm, 夜间 蒸发 强度 有 负 值 出 现 。 
2013 年 1 月 15 日 (冬季 ) 的 主要 特征 为 :平均 气温 、 
深度 10 cm 和 20 em 处 水 温 分 别 为 -5.79 9C , 
-1.94 CH -2.18*C ;与 8 月 不 同 ,10 cm 和 20 cm 
处 水 温 接近 , 且 都 高 于 气温 ;平均 湿度 .平均 净 辐 射 
和 平均 风速 分 别 为 59. 02% 、- 14. 46 W * m^ 
0.66 m s ;平均 湿度 高 于 8 月 ,平均 净 辐 射 和 风 
速 都 低 于 8 月 ,而 且 平均 净 辐 射 为 负 值 ,表明 湖水 向 
大 气 散热 ;平均 气压 .平均 湖底 压强 分 别 为 888. 90 
hPa 和 922. 52 hPa ,都 高 于 8  ;z& AIL ER YE ZR ie 
Jy 0.5 mm, 2 A vri BE «b, fa (B Hh 3,2 个 时 期 夜间 
的 蒸发 量 相近 ,而 白天 夏季 的 蒸发 量 明 显 大 于 冬季 。 
苏 木 巴 润 吉林 气象 站 数据 采集 器 每 30 min o 


集 1 组 气象 数据 , 则 1 h 内 记录 3 组 数据 。 对 3 组 
数据 求 平均 ,可 以 计算 出 逐 小 时 尺度 的 平均 气温 \ 平 
均 相 对 湿度 平均 风速 .平均 气压 平均 净 辐 射 等 要 
素 。 降 水 量 和 鞘 发 量 则 根据 该 时 段 内 数据 变化 进行 
累加 处 理 , 得 到 总 降水 量 和 荧 发 量 。 同 理 可 以 计算 
出 空气 要 素 和 净 辐 射 的 逐日 和 逐 月 平均 值 ,以 及 降 
水 量 和 蒸发 量 的 逐日 和 逐 月 累计 值 。 湖 底 压 强 数据 
每 2h 采 集 1 次 ,也 可 以 计算 出 逐日 和 逐 月 时 间 尺 
度 的 平均 值 ,进行 气压 校正 可 以 得 到 对 应 时 间 尺 度 
的 湖水 位 变化 。 考 虑 到 E601 型 蒸发 贞 置 于 湖 中 、 
蒸发 环境 最 接近 于 湖面 环境 ,本 文 将 蒸发 亚 的 蒸发 
量 作为 湖水 莹 发 量 的 近似 测量 结果 。 


3 ” 逐 小 时 蒸发 量 与 气象 要 素 的 关系 


3.1 两 两 相关 分 析 

逐 小 时 蒸发 量 与 各 种 气象 要 素 的 逐 小 时 平均 值 
相关 性 如 图 3 所 示 。 根 据 相关 系数 ,蒸发 与 净 辐 里 
的 相关 性 最 好 (R =0. 053 ) ; 与 气温 ,水 温和 风速 的 
关系 次 之 ,相关 系数 分 别 为 0. 020 .0. 022 和 0. 020; 
与 气压 ,湿度 和 风向 的 相关 性 较 弱 (RR <0.01)。 总 
体 而 言 , 逐 小 时 蒸发 量 与 各 种 同期 气象 要 素 的 相关 
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Fig.3  Correlations between hourly evaporation and meteorological factors 


性 都 较 弱 。 
3.2 多 元 线性 回归 模型 

将 气温 、 深 度 10 cm 和 20 cm 处 水 温 、 湿 度 、 气 
压 、 风速 . 净 辐 射 和 蒸发 量 数据 进行 归 一 化 处 理 , 使 
用 SPSS 统计 分 析 工 具 软 件 建立 蒸发 量 与 气象 要 素 
的 多 元 线性 回归 模型 。 本 文 参考 前 人 研究 方法 ， 
采用 逐步 回归 分 析 法 进行 变量 引入 ,根据 方差 分 析 


结果 选择 符合 判 据 条 件 且 对 因 变 量 影响 最 大 的 自 变 
量 进入 回归 模型 。 根 据 向 前 选择 法 则 进入 自 变量 ; 


然后 根据 向 后 剔除 法 ,将 模型 路 检验 概率 最 小 的 
且 符 合 剔除 判 据 条 件 的 自 变 量 剔除 ,重复 进行 直到 
模型 中 的 自 变量 均 符 合 进入 的 判 据 条 件 ,模型 外 的 
自 变量 都 不 符合 进入 模型 的 判 据 为 止 。 逐 步 回 归 分 
析 法 会 从 相关 性 很 强 的 自 变 量 中 选择 最 显著 的 因子 
放 入 回归 模型 ,避免 多 重 共 线 性 问题 。 本 文采 用 玉 
检验 的 概率 作为 判断 标准 ,使 用 系统 默认 进入 概率 
为 小 于 0. 05 ,移出 概率 为 大 于 0.1。 结果 表明 模型 
的 决定 系数 RR =0.07 ,说 明 模 型 拟 合 度 较 差 。 以 显 
著 性 水 平 小 于 0.05 作为 自 变 量 ,对 因 变 量具 有 显著 
影响 的 标准 ,得 到 多 元 线性 回归 模型 为 : 
E, = -0.047 x T, +0. 076 x T, 4 +0. 118 x R, + 
0. 073 x W, +0. 382 (1) 


IP LES 表示 逐 小 时 莹 发 强度 (mm) , T, o PORE 20 
em 处 的 湖水 温度 (% ) 。 模 拟 效果 采用 Nash. 效率 


系数 来 评价 ”“ 。 在 逐 小 时 尺度 上 ,10 cm 深度 水 
温 T; wo 与 20 em 深度 水 温 T, 之 间 的 相关 系数 很 大 
(超过 0.8) ,但 Tsa 的 作用 更 显著 ,因此 模型 中 剔除 
T T; yo 

2t AUCH ZR, NIST ZR A2 i YE Cil 20 em 
处 水 温 \ 净 辐射 和 风速 的 线性 回归 值 。 将 归 一 化 数 
据 还 原 为 真实 值 , 利 用 回归 模型 得 到 的 蒸发 量 模拟 
值 与 蒸发 由 蒸发 量 观测 值 比较 ,Nash 效率 系数 仅 为 
NSE 20.07 ,效果 较 差 。 这 表明 逐 小 时 莱 发 量 并 不 
单纯 依赖 于 同期 的 气象 要 素 , 可 能 与 更 长 时 间 尺 度 
的 气象 环境 变化 有 关 。 


4 逐日 尺度 动态 分 析 


4.1 逐日 蒸发 量 与 气象 要 素 的 相关 性 

逐日 蒸发 量 与 同期 气象 要 素 的 相关 性 如 图 4 所 
示 。 蒜 发 量 与 各 气象 要 素 的 两 两 相关 系数 普遍 大 于 
0.1 ,但 都 小 于 0.5。 其 中 ZA Et 5 ^l Kil B9 
关系 数 都 接近 0.4, 
4.2 逐日 蒸发 量 的 回归 模型 

与 3.2 所 用 方法 类 似 , 采 用 逐步 回归 分 析 法 对 
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归 一 化 的 逐日 蒸发 量 与 气象 数据 建立 多 元 线性 回归 
模型 ,得 到 的 结果 表示 为 : 

E, 20. 296 x T, 4, -0. 09 x Ry +0. 196 x 

R, +0. 297 x W, -0.078 (2) 

该 式 说 明 ,在 逐日 尺度 下 ,蒸发 量 与 20 em 处 水 温 、 
湿度 \ 净 辐射 和 风速 关系 显著 。 气 温 、7; ,6 和 7 之 
间 具 有 很 高 的 互相 关系 数 ,为 了 避免 多 重 共 线 性 问 
题 ,只 有 7 被 保留 在 回归 模型 中 。 

将 归 一 化 数据 还 原 , 利 用 回归 模型 得 到 的 蒸发 
量 模拟 值 与 蒸发 自 蒸 发 量 观测 值 比较 ,结果 如 图 5。 


明 回 归 模 型 总 体 效果 较 好 。 不 过 , 当 蒸 发 量 较 大 时 
(> 8 mm) ,线性 回归 模型 的 结果 总 体 偏 小 ,可 能 意 
味 着 蒸发 量 对 气象 要 素 的 依赖 性 是 非 线性 的 。 
4.3 非 线性 空气 动力 学 模型 
在 研究 开放 水 体 蒸 发 时 ,Mecjannet 4&0? 提出 过 

一 个 考虑 水 体面 积 的 风速 函数 ,结合 气压 差 得 到 空 
气动 力学 模型 来 估算 逐日 莱 发 量 , 即 

E,2(2.36 41.67 xW,) xA? x(e —e,) (3) 


式 中 :4 为 开放 水 体面 积 (对 本 文 观测 的 湖泊 有 4 = 
1.24 x 10* m^) ; e, 为 水 面 饱 和 水 汽 压 (kPa);e, 为 


相关 系数 R^ =0. 59, Nash 效率 系数 NSE =0.59 ,说 2 m 处 的 实际 空气 水 汽 压 (kPa) 。 直 接 使 用 式 (3 ) 进 
s [ y=0.152x+1.869 e is | y=0.184x+1.138 。 ， o i [y20.187x+1.315 P 
* R-0400 。 $ ai o | ?=0. pe- 


R?=0.399 . A E A 


E 
B 5 
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W -5 
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Fig.4  Correlations between daily evaporation and meteorological factors 
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Fig. 5 Dynamic variation of daily evaporation (a) and comparison (b) between the results estimated by 


model and the observed results 
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Fig.6 | Comparison between the daily evaporation results 
estimated by nonlinear aerodynamic model and 


the observed results 


行 模拟 ,相关 系数 达到 R = 0. 59 ,但 Nash 效率 系数 
只 有 NSE =0.26。 效 率 低 的 原因 ,可 能 主要 是 风速 
函数 与 本 地 情况 不 适应 造成 的 。 为 此 ,本 文 尝试 将 
上 述 公式 修正 为 : 

E,=(4.56+1.82 xW,) xA? x(e,-e,) (4) 


式 (4) 可 以 较 好 地 模拟 实测 蒸发 量 (图 6) , Nash 效 
率 系 数 显 车 提高 到 NSE = 0.58。 当 蒸发 量 较 大 
(>8 mm) 时 ,模拟 结果 偏 小 的 情况 得 到 改善 。 
4.4 使 用 Penman 公式 的 效果 

在 气象 水 文学 领域 ,使 用 Penman 公式 计算 逐 
日 蒸发 量 是 普遍 做 法 。 考 虑 到 这 一 点 ,本 文 将 观测 
的 逐日 蒸发 量 与 Penman 公式 "1 计算 结果 作对 比 ， 
得 到 图 7。 两 者 相关 系数 RR =0. 41, Nash 效率 系数 
NSE = 2( 负 值 表 示 模 型 失效 ) Penman 公式 计算 
的 逐日 葵 发 量 夏季 波动 大 冬季 波动 小 ,而 观测 值 全 
年 都 表现 出 较 强 的 振荡 。 而 且 ,Penman 公式 在 夏季 


一 一 观测 值 
一 一 Penman 公 式 


燕 发 量 /mm 
D o 


e 


2012-07-10 2012-12-07 2013-05-06 2013-10-03 2014-03-02 2014-07-30 


时 间 


计算 结果 偏 大 ,而 在 冬季 计算 结果 偏 低 。Penman 公 
式 给 出 的 日 平均 莹 发 量 达 到 5. 97 mm, 显著 大 于 观 
测 的 日 平均 莹 发 量 (3.74 mm) 。 莹 发 耿 观 测 值 的 确 
可 能 是 偏 低 的 ,因为 存在 一 定 的 盐 度 效应 。 蒸 发 严 
的 桶 深 是 60 cm, 两 次 换 水 期 间 的 累计 蒸发 深度 小 
于 20 cem。 根 据 前 人 文献 ,2011 年 以 来 苏 木 巴 润 
吉林 湖水 总 溶解 性 固体 (TDS) 季节 性 变化 可 能 区 间 
为 110 ~160 g - L^' , di IE dz nt 28 Ac DL p KIE n 
TDS 增加 幅度 最 大 可 达 50 g L, WEP sg sd 
观测 研究 发 现 , 在 夏秋 季 TDS 每 增加 100 g+ L7, 
可 导致 月 落 发 量 降低 49 ~ 82 mm。 若 按 比例 来 评 
佑 , 则 盐 度 效应 导致 本 文 燕 发 量 观测 的 误差 在 1.5 
mm * d 以内。 这 说 明 即 使 考虑 盐 度 效应 , 燕 发 思 
逐日 观测 值 偏 低 的 幅度 也 不 会 大 于 1.5 mm ,解释 不 
了 Penman 公式 结果 与 观测 值 如 此 大 的 偏差 。 尽 管 
在 模拟 天 然 水 体 蒸发 量 的 方法 中 ,Penman 公式 的 应 
用 最 为 广泛 ,但 笔者 推测 ,该 公式 可 能 并 不 适用 于 巴 
丹 吉林 沙漠 湖泊 的 特殊 情况 。 
4.5 逐日 蒸发 量 与 湖水 位 变化 的 相关 分 析 

从 机 理 上 讲 ,蒸发 是 驱动 湖水 位 下 降 的 一 个 因 
素 。 为 定量 评价 两 者 的 关系 ,本 文 绘制 了 逐日 燕 发 
量 与 湖水 位 变化 (上 升 为 正 、 下 降 为 负 ) 的 相关 图 
(图 8) 。 两 者 的 相关 系数 为 0.027 ,总体 相关 性 较 
35 ,但 负 相 关 的 关系 符合 湖水 位 变化 的 控制 机 理 。 

上 述 湖水 位 逐日 变化 与 逐日 蒸发 量 的 关系 可 以 
近似 表示 为 : 

E, = -0.044 x AH +3. 623 (5) 

式 中 :4 五 为 湖水 位 变化 量 。 该 模型 的 Nash 效率 系 
数 NSE =0.027 ,效果 不 佳 。 这 说 明 在 逐日 尺度 上 ， 
莹 发 量 不 能 简单 的 由 当日 湖水 位 变化 来 反映 ,还 需 


E5/mm 


图 7 观测 期 逐日 莹 发 量 实测 值 和 Penman 公式 结果 的 动态 曲线 (a) 及 数据 对 比 (b) 


Fig.7 Dynamic variation of daily evaporation (a) and comparison (b) between the results estimated by Penman Equation and 


the observed results 
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图 8 逐日 莹 发 量 与 湖水 位 变化 的 相关 性 
Fig.8 Correlation between the daily evaporation observed by 


evaporation pan and the variation of lake water level 
要 考虑 其 他 动力 学 因素 。 

` = + e E 
5 逐 月 动态 分 析 与 年 蒸发 量 


5.1 逐 月 蒸发 量 与 气象 要 素 的 相关 性 

在 月 时 间 尺 度 下 ,蒸发 量 与 气象 要 素 的 相关 性 
如 图 9 所 示 。 燕 发 量 与 气温 和 水 温 的 相关 性 最 好 ， 
相关 系数 达到 0.9 以 上 ;与 气压 \ 净 辐射 和 风速 的 相 
关 性 次 之 ,相关 系数 在 0.59 ~0.72; 与 湿度 、 风 向 相 
关 性 较 差 ,相关 系数 小 于 0.5, 
5.2 非 线 性 空气 动力 学 模型 

由 于 数据 量 较 少 ( 数 据点 不 足 20 个 ) ,不 宜 进 
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行 气象 要 素 与 逐 月 蒸发 量 的 多 元 线性 回归 分 析 。 不 
过 ,如 式 (3) 和 式 (4) 所 示 的 非 线性 空气 动力 学 模型 
所 需要 的 气象 要 素 较 少 ,利用 有 限 的 数据 校正 风速 
函数 ,可 以 得 到 : 

E, =(74.04+7.34 xW,) xA 9 x (e, -e,) (6) 


该 模型 结果 的 相关 系数 R = 0. 91, Nash 效率 系数 
NSE =0. 90 ,能 够 较 好 地 模拟 逐 月 蒸发 量 的 变化 。 
5.3 ”关于 年 燕 发 量 的 讨论 

对 上 述 逐 月 蒸发 量 进行 昧 加 ,可 以 得 到 逐年 蒸 
发 量 。 在 2013 年 6 月 至 2014 年 5 月 这 个 观测 期 
内 ,个 别 缺 失 莹 发 数据 ,根据 相 邻 的 数值 进行 线性 插 
值 。 依 据 蒸 发 亚 观 测 结 果 和 非 线性 空气 动力 学 模 
型 ,计算 得 到 的 年 莹 发 量 分 别 为 1 260 mm 和 1 312 
mm, 远 小 于 前 人 研究 ”中 推测 的 2 000 mm 以 上 
的 年 蒸发 量 。 由 于 存在 盐 度 效应 ,相对 于 实际 的 湖 
EREA RLRE EA ho BER 
也 每 次 换 水 之 前 ,TDS 的 增加 值 是 从 零 线性 增加 到 
50 g - L^' , tj EU 77 初步 可 以 用 TDS 高 于 湖 
水 25 g- L 的 程度 来 评估 , 则 夏秋 季 6 个 月 产生 的 
累计 蒸发 量 误差 小 于 120 mm, 冬 春季 蔡 发 强度 很 
小 ,附加 的 误差 应 小 于 100 mm, 因 此 , 落 发 见得 到 的 
年 蒸发 量 可 能 偏 低 100 ~ 200 mm ,按照 这 个 估计 ， 
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Fig.9 Correlation between the monthly evaporation observed by evaporation pan and the meteorological factors 
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则 实际 年 蒸发 量 很 可 能 是 1 300 ~1 500 mm ,也 显著 
低 于 2 000 mm。 沈 士 平 呈 利用 2012 年 巴 丹 吉林 沙 
漠 蒸 发 阵 观 测试 验 数据 分 析 得 到 的 全 年 淡水 营 发 量 
是 1 497 mm ,本文 对 年 落 发 量 的 评估 结果 与 此 很 接 
近 , 说 明 低 于 1500 mm 的 年 蒸发 量 对 于 巴 丹 吉林 沙 
漠 的 湖泊 来 说 是 更 加 符合 实际 的 。 


6 结论 


(1) 通过 对 沙漠 腹地 湖泊 监测 的 长 时 间 序 列 气 
象 数据 的 分 析 ,得 出 由 逐 小 时 .日 到 月 尺度 ,蒸发 量 与 
气象 要 素 的 相关 性 逐渐 增强 ,说 明 随 着 时 间 尺 度 的 延 
长 ,一些 影响 蒸发 量 的 因素 其 作用 逐渐 被 抵消 ;在 逐 
小 时 尺度 下 , 葵 发 量 与 气象 要 素 的 相关 性 都 较 低 ,日 
和 月 尺度 下 , 莹 发 量 与 气温 和 水 温 的 相关 性 最 高 。 

(2) 在 逐 小 时 尺度 下 ,多 元 线性 回归 模型 和 空 
气动 力学 模型 都 不 能 较 好 地 模拟 湖泊 水 面 蒸 发 量 ， 
表明 短 时 间 尺 度 蒸发 机 理 更 为 复杂 , 需 深 入 研究 。 
逐日 尺度 下 ,多 元 线性 回归 模型 可 以 较 好 地 模拟 湖 
itk it ze Ak, Horn KE 20. em 处 水 温和 净 辐 射 在 
蒸发 中 所 占 的 权重 较 大 ;依据 湖水 位 逐日 降幅 不 能 
很 好 地 模拟 蒸发 量 ,还 需要 考虑 其 他 因素 ;由 于 其 地 
区 性 特点 ,原始 Penman 公式 模拟 的 湖泊 水 面 蒸 发 
量 效果 较 差 。 逐 日 和 逐 月 尺度 下 ,修改 的 非 线性 空 
气动 力学 模型 可 以 都 较 好 地 模拟 湖泊 水 面 蒸发 量 。 

(3) 依据 2013 年 6 H Æ 2014 年 5 月 蒸发 血 观 
测 结果 ,并 进行 盐 度 效应 校正 ,得 到 湖泊 水 面 的 年 蒸 
发 量 为 1 300 ~1 500 mm, fHZ& Ill P3 B5 7K Hd 
量 与 湖泊 水 体 热 通 量 对 蒸发 的 影响 存在 差异 ,增加 
了 蒸发 严 测 量 薰 发 的 不 确定 性 ,值得 今后 深入 研究 。 
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Dynamic Relationship between Lake Surface Evaporation and 
Meteorological Factors in the Badain Jaran Desert 
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Abstract: 


vaporation is a key factor in lacustrine water cycle and even an entire desert. In this paper ,the dynamic relationships 


A large number of lakes are distributed in the Badain Jaran Desert located in northwest China. Water e- 


between evaporation and meteorological factors were analyzed based on the long-term meteorological and evaporative 
data monitored by automatic weather station at the Sumu Barun Jaran Lake from September 11,2012 to March 26, 
2013 and from May 8,2013 to June 23,2015 as well as the pressure intensity data at the lake bottom observed from 
September 13,2012 to June 20,2014. The multiple linear regression model , nonlinear aerodynamic model and Pen- 
man equation were used to simulate and analyze the hourly,daily and monthly evaporations. The results indicated 
that the correlations between evaporation and meteorological factors were enhanced with the time prolonging. At the 
hourly time scale,the correlations between evaporation and meteorological factors were poor,and the evaporation 
could not be simulated well by the multiple linear regression model. At the daily and monthly time scales , evapora- 
tion was closely correlated with air temperature and water temperature in the same period , and the evaporation could 
be modeled well by the multiple linear regression model and nonlinear aerodynamic model. The daily evaporation 
modeled by Penman equation was higher in summer but lower in winter. The results could provide a scientific basis 
for assessing the water consumption by lake surface evaporation in desert. 
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